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Abstract. The cause for the observed acceleration in the expansion of the universe is unknown, and 
referred to as "dark energy" for convenience. Dark energy could be an unknown energy component, 
or a modification of Einstein's general relativity. This dictates the measurements that are optimal in 
unveiling the nature of dark energy: the cosmic expansion history, and the growth history of cosmic 
large scale structure. I will examine Type la supernovae and galaxy clustering as dark energy probes, 
and discuss recent results and future prospects. 
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1. INTRODUCTION 

Solving the mystery of the observed cosmic acceleration [1] is one of the most important 
challenges in cosmology today. Current observational data are not sufficient for differen- 
tiating two likely explanations for the observed cosmic acceleration: an unknown energy 
component (dark energy, e.g., [2]), and the modification of general relativity (modified 
gravity, e.g., [3, 4]). Refs.[5, 6] contains reviews with more complete lists of references 
of theoretical models. 

The evidence for cosmic acceleration has strengthened over time. The expansion 
history of the universe is described by the Hubble parameter, H(t) = (dlna/dt) = 
d/a, where a(t) is the cosmic scale factor, and t is cosmic time. The cosmological 
redshift, z = l/a{t) — 1, is usually used as the indicator for cosmic time, because it 
can be measured for a given astrophysical object. Fig.l (left panel) shows the Hubble 
parameter H(z), as well as a, measured from current observational data [7]. Fig.l (right 
panel) shows the corresponding dark energy density function px(z) [7]. A cosmological 
constant (and no modification of gravity) continues to be consistent with data, but the 
uncertainties are large (see Fig.l, right panel). Given our ignorance of the true nature 
of dark energy, and the theoretical difficulties of explaining a tiny but non-vanishing 
cosmological constant using known physics, we need to be open minded in exploring 
and constraining alternative explanations. 

Type la supernovae (SNe la), galaxy clustering (GC), and weak lensing (WL) are 
generally considered the most powerful observational probes of dark energy. SNe la 
provide a measurement of H(z) derived from luminosity distances of SNe la. GC 
provides direct measurements of H(z), Da(z), as well as the growth rate of cosmic large 
scale structure, f g (z). WL provides measurements of H(z) and the growth factor G(z) 
(related to f g (z) via a derivative) [12]. It is important to utilize all three methods, as they 
have different systematic errors. Furthermore, GC and WL constrain different aspects in 
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FIGURE 1. Left panel: Expansion history of the universe measured from current data [7]. Data used: 
Cosmic microwave background anisotropy (CMB) data from WMAP 7 year observations [8]; 472 Type 
la supernovae (SNe la) (compiled by [9], including data from the Supernova Legacy Survey (SNLS) [9], 
the Hubble Space Telescope (HST) [10], as well as nearby SNe la [36]); galaxy clustering measurements 
from Sloan Digital Sky Survey (SDSS) Luminous Red Galaxies (LRGs) [11], 69 Gamma Ray Bursts [14], 
and the latest Hubble constant (Hq) measurement using the Hubble Space Telescope (HST) [16]. Right 
panel: Dark energy density as a function of redshift measured from the same data. 

the modification of gravity; both are required to achieve a robust test of gravity. 

Clusters of galaxies provide an independent and complementary method to probe dark 
energy [13]. Other methods, e.g., using gamma-ray bursts, old red galaxies, or radio 
galaxies [15] provide additional cross-checks on dark energy constraints. CMB data and 
independent measurements of Hq are required to break the degeneracy between dark 
energy and cosmological parameters (see e.g. [17, 7]), hence are important as well in 
constraining dark energy. 

Due to the page limit, I will focus on SNe la and GC, as these have yielded the 
strongest direct constraints on dark energy to date. I will first discuss the simplest 
general guidelines for probing dark energy, then SNe la and GC as dark energy probes 
respectively, and conclude with a brief summary of current status and future prospects 
of dark energy observational projects. 

2. PROBING DARK ENERGY AND TESTING GRAVITY 

Because of the existence of two possible explanations, dark energy and modified gravity, 
for the observed cosmic acceleration, it is critical for us to recognize that we need to 
measure two functions of cosmic time from observational data: the expansion history of 



the universe, H(z), and the growth rate of cosmic large scale structure, f g (z). Modified 
gravity models can give identical H(z) as a dark energy model by design, but the growth 
rate f g (z) is likely different in these models compared to dark energy models. The precise 
and accurate measurement of H(z) and f g (z) from observational data will allow us to 
probe the true nature of cosmic acceleration [19]. 

Dark energy is often parameterized by a linear equation of state wx (a) = wq + w a ( 1 — 
a) [20]. Because of our ignorance of the nature of dark energy, it is important to 
make model-independent constraints by measuring the dark energy density Px(z) [or 
the expansion history H(z)] as a free function of cosmic time. Measuring px{z) has 
advantages over measuring dark energy equation of state wx(z) as a free function; px(z) 
is more closely related to observables, hence is more tightly constrained for the same 
number of redshift bins used [21, 22]. Note that px(z) is related to wx(z) as follows 
[21]: 

^-expi/WHl+^l (1) 



Px(0) ~ r U<) 

Hence parametrizing dark energy with wx(z) implicitly assumes that px(z) does not 
change sign in cosmic time. This precludes whole classes of dark energy models in 
which px(z) becomes negative in the future ("Big Crunch" models, see [23] for an 
example)[24]. If the present cosmic acceleration is caused by dark energy, 

E(z) = ^= [Cl m (l+z) 3 + £l k (l+z) 2 + £l x X(z)} l/ \ (2) 

where X(z) = px(z)/px(0). Hq = H(z = 0) is the Hubble constant. Q. m and Clx are the 
ratios of the matter and dark energy density to the critical density p® = 3Hq/ '(SnG), 
and Qyt = —k/HQ with k denoting the curvature constant. Consistency of Eq.(2) at 
z = requires that Q. m + Q. k + Q. x = 1. Once E(z) is specified, the evolution of matter 
density perturbations on large scales, 8^ l \x,t) = D\(t)8(x) is determined by solving 
the following equation for D\ = 8^ (x,t) / 8(x), 1 

D'1(t)+2E(z)D\(t) - 3 -£i m {\ +z) 3 D l = 0, (3) 

where primes denote d/d(Hot). The linear growth rate 

f g (z) = dlnDi/dlna. (4) 

In the simplest alternatives to dark energy, the present cosmic acceleration is caused 
by a modification to general relativity. Ref.[26] contains examples of studies of obser- 
vational signatures of modified gravity models. A worked example is the DGP gravity 



Note that we have assumed that dark energy and dark matter are separate, which is true for the vast 
majority of dark energy models that have been studied in the literature. If dark energy and dark matter 
are coupled (a more complicated possibility), or if dark energy and dark matter are unified (unified dark 
matter models), Eq.(3) would need to be modified accordingly. Ref.[25] found the first strong evidence for 
the separation of dark energy and dark matter by ruling out a broad class of so-called unified dark matter 
models. They showed that these models produce oscillations or exponential blowup of the dark matter 
power spectrum inconsistent with observation. 



model [4, 27], which can be described by a modified Friedmann equation 2 . The right 
panel of Fig. 13 shows a DGP model that gives identical H(z) as a dark energy model, 
but gives significantly different f g (z) [19]. 

3. TYPE IA SUPERNOVAE AS DARK ENERGY PROBE 

The use of Type la supernovae (SNe la) is the best established method for probing dark 
energy, since this is the method through which cosmic acceleration has been discovered 
[1]. This method is independent of the clustering of matter 3 , and can provide a robust 
measurement of H(z) [30] through the measured luminosity distance as a function of 
redshift, cIl(z) = (1 +z)r(z), where the comoving distance r{z) from the observer to 
redshift z is given by 

r(z)=cH 1 \a k \~ l / 2 sinn[\a k \ l / 2 r(z)], (5) 
m = [^ry E(z)=H(z)/H 

where sinn(x) = sin(x), x, sinh(x) for Q. k <0,£l k = 0, and Q. k > respectively. 

Type la as standard candles. A SN la is a thermonuclear explosion that completely 
destroys a carbon/oxygen white dwarf near the Chandrasekher limit of 1.4 M & . This is 
the reason SNe la are so uniform in peak luminosity. The first challenge to overcome 
when using SNe la as cosmological standard candles is properly incorporating the 
intrinsic scatter in SN la peak luminosity. The usual calibration of SNe la reduces 
the intrinsic scatter in SN la peak luminosity (Hubble diagram dispersion) to about 
0.16 mag [31, 32]. The calibration techniques used so far are based on one observable 
parameter, the lightcurve width, which can be parametrized either as Am^ (decline 
in magnitudes for a SN la in the first 15 days after 5-band maximum, see [31]), or a 
stretch factor (which linearly scales the time axis, see [33]). The lightcurve width is 
associated with the amount of 56 Ni produced in the SN la explosion, which in turn 
depends on when the carbon burning makes the transition from turbulent deflagration 
to a supersonic detonation [34]. There may be additional physical parameters asso- 
ciated with SN la lightcurves or spectra [35] that can further improve the calibration 
of SNe la. Fig. 2 (left panel) shows a historic example of the homogeneity of SNe la [36] . 

Systematic effects of SNe la as dark energy probe. The main systematic effects for 
using SNe la to probe dark energy are: extinction by normal [37] or gray dust [38] 4 , 



2 The validity of the DGP model has been studied by [28]. 

3 Galaxy peculiar velocities from large-scale supernova surveys can be used to probe dark energy [29]. 

4 Gray dust, consisting of large dust grains, is difficult to detect by its reddening and could mimic the effect 
of dark energy [38]. Gray dust can be constrained quantitatively by the Cosmic Far Infrared Background 
[40], with no evidence found in favor of gray dust so far. Supernova flux correlation measurements can 
be used in combination with other lensing data to infer the level of dust extinction, and provide a viable 
method to eliminate possible gray dust contamination in SN la data [41]. 
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FIGURE 2. Left panel: Hubble diagrams showing 26 SNe la with B max — V max < 0.20 from the 
Calan/Tololo sample [36]. This sample provided half of the data for the discovery of the cosmic ac- 
celeration in 1998 [1]. The solid lines indicate Hubble's law; perfect standard candles (with a = 0) fall on 
these lines. Right panel: Hubble diagrams of SNe la in the NIR bands. Note that these SNe la have only 
been corrected for dust extinction; no corrections have been made for lightcurve width. [46] 



weak lensing amplification by cosmic large scale structure [39], and possible evolution 
in the peak luminosity of SNe la. 

Recent data show that the apparent dust extinction of SNe la is very different from 
the typical extinction law due to Milky Way dust, possibly due to the mixing of intrinsic 
SN la color variation with dust extinction, or variations in the properties of dust [42]. 
The extinction by dust can be corrected using multi-band imaging data, especially near 
infrared (NIR) observations of SNe, since dust extinction decreases with wavelength. 

The weak lensing amplification of SNe la by cosmic large scale structure can be 
modeled by a universal probability distribution function for weak-lensing amplification 
based on the measured matter power spectrum [43]. The effect of weak lensing on the 
SN la data can be minimized through flux- averaging [44]. Figs. 3 shows the 2D marginal- 
ized contours of (w,Q m ,^#) (where is a nuisance parameter), assuming a constant 
equation of state for dark energy, w, and a flat universe. Note that the inclusion of sys- 
tematic errors of SNe leads to significantly larger uncertainties in estimated parameters, 
compared to when only statistical errors of SNe are included [7]. Clearly, flux- averaging 
(thick solid lines) leads to larger errors on dark energy and cosmological parameters if 
only SN la data are used. However, when other data are added, flux- averaging leads to 
smaller errors on dark energy (see Figs. 4-5) because flux- averaging increases the con- 
cordance of SNe la with other data. 

The evolution in SN la peak luminosity could arise due to progenitor population drift, 
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FIGURE 3. The 2D marginalized contours of (w,Q. m ,^) for SNe data compiled by [9] (with and 
without flux-averaging), assuming a fiat universe. The contours are at 68% and 95% confidence levels. [7] 
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FIGURE 4. The 2D marginalized contours of (wo,w a ) and (wo, wo. 5) for SNe data (with and without 
flux-averaging) combined with galaxy clustering (CW2), CMB, Hq, and GRB data (same data as in Fig.l). 
The contours are at 68% and 95% confidence levels. [7] 
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FIGURE 5. The 2D marginalized contours of (Xo ^,Xo£7,Xi.Q,Cl m ,Q.k) f° r SNe data (with and without 
flux-averaging) combined with galaxy clustering (CW2), CMB, Ho, and GRB data (same data as in Fig.l). 
The contours are at 68% and 95% confidence levels. [7] 



0.5 



0.0 




V R I 

Wavelength Bond 



FIGURE 6. Dispersion in peak magnitude (measured at the first lightcurve maximum) as a function of 
wavelength band for SN la models with 56 Ni masses between 0.4 and 0.9 M (Kasen 2006 [49]). 



since the most distant SNe la come from a stellar environment very different (a much 
younger universe) than that of the nearby SNe la. However, with sufficient statistics, 
we can subtype SNe la and compare SNe la at high redshift and low redshift that are 
similar in both lightcurves and spectra, thus overcoming the possible systematic effect 
due to progenitor population drift [45]. 

Optimized observations of SNe la. NIR observations of SNe la provide additional 
strong advantages beyond being relatively dust-free. SNe la are better standard candles 
at NIR wavelengths compared to the optical wavelengths [46, 47, 48]. The right panel of 
Fig. 2 shows the Hubble diagram of SNe la in the NIR, without the usual lightcurve width 
correction. The smaller intrinsic dispersion of SN la peak luminosity in the NIR can be 
explained by the theoretical modeling of SN la lightcurves using time-dependent multi- 
group radiative transfer calculations (Kasen 2006 [49]). Fig. 6 shows the dispersion in 
peak magnitude (measured at the first lightcurve maximum) as a function of wavelength 
band for SN la models with 56 Ni masses between 0.4 and 0.9 M [49]. 

It is important to obtain high quality spectra of SNe la (including NIR spectra, see 
[50]), since the spectra of SNe la have been shown to provide calibration relations that 
decrease the scatter of SNe la in the Hubble diagram, and make SNe la better distance 
indicators. The correlation between SN la spectroscopic features and luminosity has 
been found in the observational data (see, e.g., [51]). More recently, Bailey et al. (2009) 
[52] used the Nearby Supernova Factory spectrophotomery of 58 SNe la to perform an 
unbiased search for flux ratios that correlate with SN la luminosity. They found that the 
642/443 nm flux ratio is most strongly correlated with SN la absolute magnitudes. The 
correlation of SN la spectroscopic features and luminosity can be understood through 
comparing theoretical modeling with observational data. Hachinger et al. (2008) [53] 
found that the strength of the Si II A 5972 line may be a very promising spectroscopic 




FIGURE 7. The comparison of an ultra deep supernova survey [55] with a much shallower survey in 
the reconstruction of the dark energy density px(z) as a free function of cosmic time [54]. 



luminosity indicator for SNe la, with the correlation between Si II A 5972 strength and 
luminosity resulting from the effect of ionization balance. 

The key to the efficient use of SNe la for probing dark energy is to obtain the largest 
possible unbiased sample of SNe la at the greatest distances from the observer [54]. This 
is achieved by an ultra deep survey of the same areas in the sky every few days over at 
least one year [55]. Given the same observational resources, an ultra deep supernova 
survey is superior to a much shallower survey. A sufficiently deep supernova survey is 
required to reconstruct the dark energy density px(z) as a free function of cosmic time 
(i.e., to measure H(z) precisely, see Fig. 7) [54]. 



4. GALAXY CLUSTERING AS DARK ENERGY PROBE 

Baryon acoustic oscillations (BAO) as a standard ruler. At the last scattering of 
CMB photons, the acoustic oscillations in the photon-baryon fluid became frozen, and 
imprinted their signatures on both the CMB (the acoustic peaks in the CMB angular 
power spectrum) and the matter distribution (the baryon acoustic oscillations in the 
galaxy power spectrum). Because baryons comprise only a small fraction of matter, 
and the matter power spectrum has evolved significantly since last scattering of photons, 
BAO are much smaller in amplitude than the CMB acoustic peaks, and are washed out 
on small scales. BAO in the observed galaxy power spectrum have the characteristic 
scale determined by the comoving sound horizon at the drag epoch (which occurred 
shortly after recombination), which is precisely measured by the CMB anisotropy data 
(see, e.g., [8]). Comparing the observed BAO scales with the expected values gives H(z) 
in the radial direction, and Da(z) = r(z)/(l +z) (the angular diameter distance) in the 
transverse direction [56, 57]. Fig. 8 shows the first detection of the BAO peak from a 
sample of the SDSS LRGs [58]. 

BAO represents only a fraction of the cosmological information contained in galaxy 
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FIGURE 8. The spherically-averaged galaxy correlation function measured from the SDSS data, clearly 
showing a peak corresponding to the BAO scale at ~ lOO/i^'Mpc [58]. 



clustering data. A flux-limited galaxy redshift survey can allow us to measure the cosmic 
expansion history H(z) through BAO only or the shape of the galaxy power spectrum 
P(k) (which includes BAO as features) measured from the galaxy distribution, and the 
growth history of cosmic large scale structure f g (z) through independent measurements 
of redshift-space distortions and the bias factor between the distribution of galaxies and 
that of matter [19]. Assuming linear bias, the combination f g (z)o% m {z) can be measured 
directly and used to test gravity [59]. 

Current GC Measurements. The first simultaneous measurements of H(z) and 
Da(z) from galaxy clustering data was made very recently by Chuang & Wang (2011) 
[11], based on a Markov Chain Monte Carlo (MCMC) analysis of the two-dimensional 
two-point correlation function (2D 2PCF) they measured from the flux-limited sample 
of LRGs from the SDSS Data Release 7 [1 1]. Fig.9 shows the 2D 2PCF measured from 
the SDSS LRGs and a single LasDamas SDSS LRG mock catalog for comparison. 
The similarity between the data and the mock in the range of scales used (indicated 
by the shaded disk) is apparent. Due to the current limitations in the modeling of 
systematic effects, only the quasi-linear scale range of s = 40 — 120/i _1 Mpc is used 
for a conservative estimate in this analysis. Chuang & Wang (201 1) [1 1] obtained 

//( z = 0.35) = 82.1+2jkms- 1 Mp C - 1 , D A (z = 0.35) = 1048+g> Mpc (6) 

without assuming a dark energy model or a flat universe. Scaling H(z) and Da(z) using 
r s(zd) (the sound horizon at the drag epoch) in the MCMC analysis, they found that the 
derived measurements of 




H(0.35)r s (z d ) = 13020 ±530km/s, r s (z d )/D A (0.35) = 0.1518 ±0.0062 (7) 
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FIGURE 9. The two-dimensional two-point correlation function (2D 2PCF) measured from SDSS 
DR7 LRGs (left panel) and a LasDamas SDSS LRG mock catalog (right panel) in a redshift range 
0.16 < z < 0.44 (solid black contours), compared to a theoretical correlation function with parameters 
close to the best fit values in the likelihood analysis (dashed red contours). In both figures, the shaded 
disk indicates the scale range considered (s = 40 — 120 /z^'Mpc ) in this study. The thick dashed blue 
circle denotes the baryon acoustic oscillation scale. The observed 2D 2PCF has been smoothed by a 
Gaussian filter with rms variance of 2/i~'Mpc for illustration in these figures only; smoothing is not used 
in our likelihood analysis. The contour levels are % = 0.5,0.1,0.025,0.01,0.005,0. The ^ = contours 
are denoted with dotted lines for clarity. [11] 



are nearly uncorrelated (with a normalized correlation coefficient of r = —0.0584), have 
tighter constraints and are more robust with respect to possible systematic effects. This 
is as expected, since //(0.35) r s (zd) and r s (zd)/DA(0.35) correspond to the preferential 
redshift separation along the line of sight, and the preferential angular separation in 
the transverse direction respectively; these in turn arise from the BAO in the radial 
and transverse directions. The measurable preferential redshift and angular separations 
should be uncorrelated since they are independent degrees of freedom. The presence of 
the BAO (although only marginally visible in Fig. 9) leads to tight and robust constraints 
on H (0.35) r s (zd) and r s (zd) /Da{0. 35). Since most of the constraining power in this 
analysis comes from fitting the overall shape of the galaxy correlation function on 
quasi-linear scales, and not from fitting the BAO peaks, these measurements are galaxy 
clustering measurements (rather than BAO only measurements). 

The constraints in Eq.(7) can be used to combine with cosmic microwave background 
and any other cosmological data sets to constrain dark energy, as no priors were imposed 
that would affect the combined constraints [60, 11]. These results have significant 
implications for future surveys in establishing the feasibility of measuring both H(z) 
and Da(z) from galaxy clustering data. 

Fig. 10 and Fig. 11 shows the first results from the WiggleZ Dark Energy Survey 
obtained by Blake et al. (201 lab) [61, 62]. Note that Fig. 10 assumes fixed values of 
(Q/,/? 2 , n$, h, Og), and the background cosmological model is assumed to be known for 
Fig. 11. When dramatically larger data sets become available from the next generation 
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FIGURE 10. Measurements of the baryon acoustic peak at redshifts z= 0.44, 0.6 and 0.73 in the galaxy 
correlation function of the final data set of the WiggleZ Dark Energy Survey. [61] 
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FIGURE 11. Measurements of /g(z)og m (z) from the WiggleZ Dark Energy Survey data. [62] 



galaxy redshift surveys, it will be possible to extract both distance and growth rate 
measurements simultaneously without imposing strong priors (see, e.g., [63]). 

Possible systematic differences in different GC measurements. Ref.[7] showed 
that there may be systematic differences in different GC measurements. Fig. 12 shows 
the 2D marginalized contours of (w, Q. mi 0.%) for different GC measurements combined 
with CMB, H , and GRB data. The first row of Fig.12 compares the H(z = 0.35)r s (z d ) 
and r s (zd) / Da(z = 0.35) measurements by Chuang & Wang (2011) [11] with their 
^0.35 = r s (zd)/Dv(z = 0.35) measurement (both from SDSS DR7 LRGs), as well as 
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FIGURE 12. The 2D marginalized contours of (w,Q.,„,Q.i ( ) for different galaxy clustering measure- 
ments combined with CMB, Ho, and GRB data. The contours are at 68% and 95% confidence levels. 



the do.2 and Jq.35 measurements by Percival et al. (2010) [64] from SDSS DR7 LRG 
and main galaxy samples and 2dFGRS, and the Jo.6 measurement by Blake et al. (201 1) 
from the WiggleZ survey [66] combined with the <io.i06 measurement by Beutler et al. 
(2011) from 6dF GRS [67]. 

For the Chuang & Wang (2011) [11] GC measurements (CW2 and CW1), the con- 
straints on w are tightened significantly by going from spherically-averaged data (CW1), 
i.e., d 035 , to 2D data (CW2), i.e., H(z = 0.35)r s (z d ) and r s (zd)/D A (z = 0.35), as indi- 
cated by comparing the thin solid contours (CW1) to thick solid contours (CW2) in the 
first row of Fig. 12. This is as expected, as more information from GC is included in 
CW2 compared to CW1 . Both the Percival et al. (2010) GC measurements (WP) and the 
combined WiggleZ survey and 6dF GRS measurements (CB+) favor w < — 1 (similar 
results were found by [65] using GC measurements from [64]), while the Chuang & 
Wang (201 1) [11] GC measurements favor w = — 1 . 

The second row in Fig. 12 compares the do.2 and do.35 measurements by Percival et al. 
(2010) [64] (WP2), with their measurements of Jo. 2 and do.35 separately. Clearly, most 
of the constraining power on w comes from c?o.35- While the do.2 measurement favors 
w = — 1, the Jq.35 measurement favors w < — 1. 

The measurements of rfo.35 by Chuang & Wang (201 1) [1 1] and Percival et al. (2010) 
[64] are similar in precision, but differ systematically: d™ 5 = r s (z c i)/Dy(z = 0.35) = 
0.1161 ±0.0034, while dfj 5 = r s (z d )/D v (z = 0.35) = 0.1097 ± 0.0036. The lower 
measured value of dff£ 5 implies a smaller H(z = 0.35), which in turn implies a more 
negative w. When combined with CMB, Hq, and GRB data, d^ 5 favors w = — 1, while 




FIGURE 13. Left panel: The average two-dimensional two-point correlation function (2D 2PCF) mea- 
sured from 160 LasDamas SDSS LRGfull mock catalogs (solid black contours), compared to a theoretical 
model with the input parameters of the LasDamas simulations (dashed red contours). The line types and 
contour levels are the same as in Fig. 9 [11]. Right panel: Current and expected future measurements of 
the cosmic expansion history H(z) = HqE(z) and the growth rate of cosmic large scale structure f g (z)). 
The future data correspond to a magnitude-limited NIR galaxy redshift survey covering > 10,000 square 
degrees and 0.5 < z < 2. [19] 



^0 35 f avors w < —1. Note that these two measurements used different methods to 
analyze GC data: Chuang & Wang (2011) used the galaxy correlation function, while 
Percival et al. (2010) used galaxy power spectrum. It is not surprising that they lead to 
different distance measurements from GC. 

Future Prospects. The data contours in the left panel of Fig. 13 gives a sense of the 
exquisite precision the galaxy 2D 2PCF can be measured when significantly larger data 
sets become available. Fig. 13 (left panel) [11] shows the averaged 2D 2PCF measured 
from 160 LasDamas mock catalogs compared with a theoretical model. The contour 
levels are apparent in the measured 2D 2PCF even though no smoothing is used (in 
contrast to the noisy current data, see Fig. 9); this is due to the reduction of shot noise 
achieved by averaging over 160 mock catalogs. Clearly, the 2D theoretical model used 
by [11] provides a reasonable fit to data on intermediate (and quasi-linear) scales. 

A flux-limited galaxy redshift survey can allow us to measure both//(z) and f g (z) [18, 
19]. The measurement of f g (z) can be obtained through independent measurements of 
redshift- space distortion parameter /3 = f g (z)/b [68] and the bias parameter b(z) (which 
describes how light traces mass) [18]. The parameter /3 can be measured directly from 
galaxy redshift survey data by studying the observed redshift-space correlation function 
[69, 70]. We can assume that the galaxy density perturbation 8 g is related to the matter 
density perturbation <5(x) as follows [71]: 8 g = b8(x) +^2<5 2 (x)/2. The galaxy bispec- 



trum is (8 gkl 8 gk2 8 gki ) = {litf {P g {k,)P g {k 2 ) [J(k h k 2 )/b + b 2 /b 2 ] + cyc.}8 D (k l + 
k2 +k3), where 7 is a function that depends on the shape of the triangle formed by (ki, 
k2, k3) in k space, but only depends very weakly on cosmology [72]. Ref.[72] developed 
the method for measuring b(z) from the galaxy bispectrum, which was applied by [73] 
to the 2dF data. Independent measurements of f}(z) and b(z) are very limited at present 
[69, 73]; this will change dramatically in the near future. 

The right panel of Fig. 13 shows how well a flux-limited NIR galaxy redshift survey 
covering >10,000 square degrees and 0.5 < z < 2 can constrain H(z) and f g (z), com- 
pared with current data [19]. The bottom half of the right panel in Fig. 13 shows the f g (z) 
for a modified gravity model (the DGP gravity model) with O.^ = 0.25 (solid line), as 
well as a dark energy model that gives the same H{z) for the same O.^ (dashed line). 
The cosmological constant model from the top half of the right panel in Fig. 13 is also 
shown (dotted line). Clearly, current data can not differentiate between dark energy and 
modified gravity. A very wide and deep galaxy redshift survey provides measurement 
of fg(z) accurate to a few percent; this will allow an unambiguous distinction between 
dark energy models and modified gravity models that give identical H(z) (see the solid 
and dashed lines in the bottom half of the right panel of Fig. 13). 

The systematic effects of BAO as a standard ruler are: bias between luminous matter 
and matter distributions, nonlinear effects, and redshift distortions [56]. Cosmological 
N-body simulations are required to quantify these effects [74]. Ref.[75] shows that 
nonlinear effects can be accurately taken into account. Ref.[76] shows that the BAO 
signal is boosted when bias, nonlinear effects, and redshift distortions are properly 
included in the Hubble Volume simulation. For a detailed discussion, see [6]. 

5. SUMMARY: CURRENT STATUS AND FUTURE PROSPECTS 

There are a large number of dark energy surveys that are ongoing, planned, or pro- 
posed. Ongoing projects include Carnegie Supernova Project (CSP), ESO Visible and 
Infrared Survey Telescope for Astronomy (VISTA) Surveys, Panoramic Survey Tele- 
scope & Rapid Response System (Pan-STARRS), Hobby-Eberly Telescope Dark En- 
ergy Experiment (HETDEX), and Sloan Digital Sky Survey (SDSS) III, and the Dark 
Energy Survey (DES) [77]. Selected future projects include Large Synoptic Survey Tele- 
scope (LSST) [78], and Euclid [63, 79]. Proposed future projects include the Wide-Field 
Infrared Survey Telescope (WFIRST), BigBOSS, and Square Kilometre Array (SKA) 
[80]. It is critical to remember that the challenge to solving the dark energy mystery 
will not be the statistics of the data obtained, but the tight control of systematic effects 
inherent in the data. A combination of all three most promising methods (SNe la, GC, 
and WL) should be used, each optimized by having its systematics minimized by de- 
sign [81]. It is an exciting time in cosmology. We can expect to make ground-breaking 
discoveries about the nature of dark energy within the next decade or two. 

I am grateful to Chris Blake, Daniel Eisenstein, Mario Hamuy, Dan Kasen, and Kevin 
Krisciunas for permission to use their figures in this proceeding paper. This work was 
supported in part by DOE grant DE-FG02-04ER41305. 
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